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The stationary thermionic electron emission currents from heated metals are compared against an 
analytical expression derived using a non equilibrium quantum Kappa energy distribution for the 
electrons. This later depends on the temperature decreasing parameter k(T) which can be estimated 
from the raw experimental data and characterizes the departure of the electron energy spectrum from 
the equilibrium Fermi-Dirac statistics. The calculations accurately predict the measured thermionic 
emission currents for both high and moderate temperature ranges. The Richardson-Dushman law 
governs the electron emission for large values of Kappa or equivalently, for moderate metal temper¬ 
atures. The high energy tail in the electron energy distribution function which develops at higher 
temperatures or lower Kappa parameters, increases the emission currents well over the predictions 
of the classical expression. This analysis also permits the quantitative estimation of the departure 
of the metal electrons from the equilibrium Fermi-Dirac statistics. 

PACS numbers: 05.90.+m 05.30.Fk 79.40.+Z 


I. INTRODUCTION 

The thermionic electron emission from metals at high 
temperatures is today on the basis of countless techni¬ 
cal applications. Different materials are currently used 
as hot cathodes to produce beams of negatively charged 
particles in electron guns, plasma sources or microwave 
devices such as klystrons or traveling wave tubes |l], and 
vacuum thermionic energy conversion devices . 

The stationary flow of thermionic electrons from the 
metal surface is currently calculated by the classical 
Richardson-Dushman (RD) equation. This model con¬ 
siders the thermal equilibrium between the electron gas 
and the metal lattice and makes use of the Fermi-Dirac 
distribution for the electron energy spectrum 0]. The 
RD thermionic electron current density Jrd{T) essen¬ 
tially relies on the work function Wf of the metal and 
its temperature T, disregarding other important factors, 
such as the geometry [4], [5] of the emitting surface or its 
physical state 6]. 

However, the electron energy distribution in metals 
at high temperatures frequently differs from the Fermi- 
Dirac statistics. The stationary RD electron emission 
regime represents the final step after different energy 
thermalization processes with shorter time scales 0, 7}- 
[9]. This time dependent decay towards the stationary 
thermal equilibrium have been studied by exposing the 
metal surfaces to ultrashort laser pulses. The subsequent 
evolution of the energy spectrum is later monitored to 
determine the electron energy relaxation rates 0i0. 

The average electron energy relaxes fast, within the 
femtosecond time scale by collisions between electrons. 
Therefore, the electron T e and metal T temperatures dif¬ 
fer within these short times scales. The thermalization 
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between the electron gas and the metal lattice requires of 
few tens of picoseconds due to the larger mass of phonons 
001- The RD regime is considered to provide an ade¬ 
quate description of the electron emission for long charac¬ 
teristic times when the electron gas and the metal lattice 
reach the thermal coupling T = T e (l, 0, ITlj |- Never¬ 
theless, this pure RD thermionic emission regime of elec¬ 
trons has not been observed because the energy of the 
laser pulse lies over damage thresholds of most metals 
0 . 

This physical description relies on the implicit assump¬ 
tion of the long time scale relaxation to a thermal equilib¬ 
rium where the Fermi-Dirac statistics describes the elec¬ 
tron energy spectrum 01 - This might be not always the 
case, even under an efficient energy transfer between the 
electron gas and metal lattice. For high temperatures the 
metal might remain out of equilibrium because additional 
energy exchange processes take place, such as the intense 
emission of electromagnetic radiation, the development 
strong surface thermal gradients or the nonuniform elec¬ 
tron emission. The high energy interactions between the 
electrons and phonons couple the electron gas with the 
overheated metal lattice, producing groups of fast elec¬ 
trons 0,0. Therefore, even for long time scales, the elec¬ 
tron energy spectrum might differ from the Fermi-Dirac 
energy distribution in hot metals. 

Consequently, the equilibrium Fermi-Dirac statistics 
needs to be replaced to account for these high energy 
electron groups. The usual equilibrium statistical physics 
cannot cope with these non-conventional electron energy 
distribution functions previously proposed in different 
fields [l2j]. Namely, the derivation of a generalized Planck 
radiation law_[13l|. the description of high temperature 
Fermi gases [l4f | or the out of equilibrium warm dense 
matter Q0. 

Recently, an analytical expression was derived for 
the thermionic current density J K (T), using a modified 
Kappa energy distribution for the electrons jl6| . This 
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Figure 1: The diagram of experimental setup where, Vh, is the 
voltage drop across the tungsten wire heated by the DC cur¬ 
rent, Ih- The dotted box represents the scheme of the sweep 
bias circuit which applies an amplified sawtooth waveform 
V p (t) to the heated wire and measure the emission current, 
I P {t), as the voltage drop, AV, across a precision resistor, Rd. 
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Figure 2: The electron emission currents 7 e (T) as a function 
of the bias potential V v for different metal temperatures. 


II. THE NON EQUILIBRIUM ELECTRON 
EMISSION MODEL 

The stationary thermionic electron current from met¬ 
als at high temperatures could be calculated using the 
stationary, quantum non-equilibrium Kappa statistics 
0, 


formula corresponds to the high energy, non equilibrium 
Fermi-Dirac distribution [17], [lg]. The index k of this 
distribution accounts for the high energy tail in the elec¬ 
tron energy spectrum, as well as the Kappa deformed 
statistics previously proposed ra¬ 
in this paper the measurements of the stationary 
thermionic emission currents from metals are found in 
good agreement with the theoretical predictions of Ref. 
1161 . The RD expression underestimates the experi¬ 

mental emission currents for high metal temperatures, 
whereas Jrd(T) and J K (T) agree for moderate values. 


( e - e F \ -(K+1) 

h(Te, E) = C K (T e ) [l + (2-1) 

Here k B is the Boltzmann constant and E is the kinetic 
energy of electrons of mass m e , which is evaluated with 
respect to the Fermi level cf- The other parameters are 
7 i = tF/ksTe and T K = (n - 3/2 + 71 ) T e . Eq. (12.11) is 
the approximation for high energies of the generalized 
non equilibrium quantum statistics of Refs. fl7l and jl8| . 
The normalization factor C K (T e ) in Eq. (12.11) scales to 
the metal electron density n eo as, 


As we shall see, the electron energy distribution might 
be approximated using the temperature dependent k(T) 
index of this statistics, which could be determined from 
the experimental data. This J K {T ) accounts for the con¬ 
tribution of high energy electrons in the thermionic elec¬ 
tron emission. 
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Additionally, it might be employed as a quantitative 
estimation of the the departure from the thermal equi¬ 
librium described by the Fermi-Dirac statistics which is 
recovered in the limit of large n. To our best knowledge, 
this present analysis is one of the few cases where the 
predictions of a non equilibrium quantum Kappa distri¬ 
bution are directly compared against raw experimental 
data. 


The low values of k in Eq. m account for a large 
high energy tail in the electron energy spectrum, whereas 
large indexes n recover the classical Fermi-Dirac statistics 

m 

For long time scales, well over the microsecond range 
equal temperatures T e = T are considered for the elec¬ 
tron gas and the metal lattice mm. The average elec¬ 
tron thermal energy proportional to k B T e is always much 
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lower than the work function, Wf. In these conditions 
the analytical expression for the stationary thermionic 
electron current density J K (T) deduced from Eq. (12.111 
reads (l6| . 
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For low k the thermionic electron fluxes predicted by Eq. 
(12.21) are higher than those calculated using the classical 
Richardson-Dushmann expression Jrd(T)' which is also 
recovered in the opposite limit of large k |16j . 

The contribution in J K (T ) of the high energy tail of 
the electron energy distribution function increases with 
the metal temperature. This effect is incorporated by 
means of a temperature dependent index k(T) in either 
f K (T,E) and J K {T). 


III. EXPERIMENTS 



Figure 3: Comparison between the experimental data 
I a (T)/I S (T 0 ) and the normalized thermionic current densi¬ 
ties jn{T) = Jk{T) / Jn. m (T 0 ) (labeled as J. Theor.) by Eq. 
(12.21) using the fitting k(T) of Fig. [3] The dashed line repre¬ 
sents the predictions of the normalized Richardson-Dushmann 
equation jnD{T) as a function of the metal temperature. 


The predictions of Eq. m were checked against 
the measurements of the thermionic electron currents 
I e (T,V p ) from a DC heated tungsten and in Fig. |T| is 
depicted a simplified scheme of the experimental setup. 
The loop shaped wire with a typical length of 22 mm 
and 0.08 mm in diameter was placed at the end of a elec¬ 
trically insulated ceramic shaft. The emitting wire was 
approximately located at the center of a cylindrical vac¬ 
uum chamber of 0.8 m in length and 0.4 m in diameter 
evacuated down to typical pressures below 10 -5 rnB of 
Argon. The elastic collisions between electrons and neu¬ 
tral is the dominant collisional process. According to [20] 
the cross sections are 9 — 20 x 10~ 2 ° m 2 , which gives typ¬ 
ical mean free paths between 20 and 44 m, much larger 
than our experimental arrangement. This fact excludes 
the production of additional charged particles because 
the ionizing and elastic collisions between electron and 
neutrals are negligible. 

The tungsten wire was heated up to thermionic elec¬ 
tron emission by a DC currents of/ft,~0.8 — 1.0A with 
voltages Vh —2 — 3V and electrically biased with respect 
to the grounded walls of the vacuum tank. For these thin 
tungsten filaments the temperature gradients along the 
hot wire could be neglected and energy power losses are 
mainly caused by radiation emission. Therefore, the tem¬ 
peratures were determined using a well known expression 
relating the DC heating current alon g th e wire with T, 
which is regarded essentially uniform |2l[-|23|. 

The sweep system impress an amplified time depen¬ 
dent sawtooth signal, V p (t), with a repetition pulse of 2 
kHz to one leg of the heated tungsten wire [24]]. As shows 
Fig. [T] the emitted thermionic electron current I e (T,V p ) 
is obtained from the current, I p (t), measured as the volt¬ 
age drop, Aid, across the precision resistance, Rd- 


The typical electron emission currents J e (T, V p ) are 
represented in Fig. [2] for different metal temperatures 
T as a function of the bias potential Vp of the wire. 
For negative bias potentials V p the thermionic electrons 
are emitted from the wire towards the vacuum chamber 
walls, whereas no electron current is observed for positive 
potentials. 

The lower currents in Fig. [2] were of tens of //A corre¬ 
sponding to T < 2100 K that raised up to the 0.1-0.2 rnA 
range for metal temperatures over 2200 K. For low posi¬ 
tive potentials where the electron emission becomes neg¬ 
ligible, I e (T,V p ) reaches a flat negative current for volt¬ 
ages over a threshold of about 5 volts. The thermionic 
electron emission takes place for negative bias potentials 
and large slopes in the range -5 to -20 volts correspond 
to the space charge effects around the wire. These are 
outweighed over a typical threshold bias potential of the 
wire corresponding to the knee of J e (T, V p ) where the flat 
electron emission current becomes weakly dependent on 
the bias potential. 

Then, the flat saturation currents I S {T) ~ I e (T,V p ) 
for V p < —20 volts past the knee of 7 e (T, V p ), measure 
the maximum emitted electron thermionic current. The 
moderate bias voltages V p involved in Fig. [2]exclude from 
this pure thermionic regime the thermofield or Fowler- 
Nordheim or field emission modes [25|, [26[ . 


IV. EXPERIMENTAL DATA ANALYSIS 


The theoretical predictions of Eq. (12.211 are compared 
with the measurements of 7 S (T) in Fig. [3] For low metal 
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temperatures (typically below 2000 K) both, the RD ex¬ 
pression and Eq. (12.21) give similar thermionic electron 
emission currents. However, the increment in the metal 
temperature leads to important differences. 

Additionally, the values of index k(T) of f K (T, E ) could 
be evaluated from the maximum thermionic electron 
emission currents I S (T) and are represented in Fig. [I] 
The low temperature limit corresponds in Fig. [3] to large 
values of the index k(T) where f K (T, E ) becomes similar 
to the equilibrium Fermi-Dirac statistics. The decreasing 
k(T) index characterizes the f K (T,E) for growing metal 
temperatures. 

The following procedure was used in Figs. [3] and [4] 
to avoid the need of the accurate determination of the 
electron emitting metal surface, as well as the particular 
value of the Richardson constant for the metallic samples. 

First, we obtain the parameter n m for the lowest metal 
temperatures T a ~ 2100 K in Fig. © which corresponds 
to the Richardson-Dushman electron emission regime. 
From Eq. m and the classical RD expression we obtain 
k m as the root of a nonlinear equation by setting, 

Jrd{T 0 ) 

This permits to evaluate the lower current density 
J Km (T 0 ) using Eq. (12.21) which corresponds to the lower 
measured thermionic current I S (T 0 ) of Fig. [2] The max¬ 
imum values is about k m ~ 25.6 — 26.0 in Fig. [3] for the 
three sets of experimental data. 

Next, using the maximum thermionic emission cur¬ 
rents I S (T) for the metal temperature T the ratios 
I S (T)/I S (T 0 ) are therefore, 


UT) = J K (T) = 

Is(T 0 ) J Km (T 0 ) Jk{ ’ 

Using Eq. m the values k(T) < n m for increasing 
temperatures T > T a could be again evaluated as the 
roots of a nonlinear equation. 

Finally, the temperature decreasing values of k(T) of 
Fig. ©could be approximated by k(T) = a—bT by means 
of a least-squares fitting and these empirical expressions 
are later introduced in Eq. (EH). This procedure allows 
to compare its theoretical predictions with the experi¬ 
mental ratios I S (T)/I S (T 0 ). The normalized curves j K (T) 
along with the experimental ratios I S (T)/I S (T 0 ) are rep¬ 
resented in Fig. [3] as well as, 


Jrd(T) 


Jrd(T) 

Jrd(T 0 ) 


which it were calculated by using the classical 
Richardson-Dushmann equation. 


V. CONCLUSIONS 

The classical Richardson-Dushman equation based on 
the equilibrium Fermi-Dirac statistics adequately de¬ 



Figure 4: The parameter k(T ) corresponding to the exper¬ 
imental data fitting, using the least-squares method, of the 
three data sets of Fig. © to k(T) = a — bT, where T is the 
metal temperature. 


scribes the thermionic electron emission for low and mod¬ 
erate metal temperatures. However, additional energy 
exchange processes bring the metal lattice far form the 
thermal equilibrium at higher temperatures. Thus, a 
fraction the lattice energy is also transferred to the elec¬ 
tron gas for long time scales. In these conditions, the elec¬ 
tron energy spectrum differs from the Fermi-Dirac statis¬ 
tics and could be approximated by the effective, non equi¬ 
librium quantum Kappa distribution of Eq. (12.11) . This 
later recovers the Fermi-Dirac statistics for low temper¬ 
atures or equivalently, for large values of the «(T) index 
in UT,E). 

The empirical determination of the index k(T) from 
the stationary thermionic electron emission currents 
J S (T) for different metal temperatures was discussed in 
Sec. IV. The maximum value K m for the lower tem¬ 
perature in Fig. U corresponds to the RD regime and 
the decreasing function k(T ) quantitatively characterizes 
the departure of the electron energy spectrum from the 
Fermi-Dirac statistics [ij|. This transition rate is gov¬ 
erned by the slope b in the empirical fit to to k(T) = 
a — bT of Fig. 3] which also characterize the increasing 
high energy tail in f K (T, E) with the metal temperature. 

Fig. [3] shows the good agreement between the exper¬ 
imental data and the thermionic electron emissions cur¬ 
rents J k (T) calculated using f K (T,E) for the electron 
energy spectrum and the empirical expressions for k(T). 
Additionally, the measurements of Fig. [3] show the in¬ 
crease of the thermionic current by orders of magnitude 
over RD law predictions for high metal temperatures. 
These increments are caused by the departure for the 
thermal equilibrium of both the metal lattice and the 
electron gas not contemplated by the equilibrium Fermi- 
Dirac statistics employed in the derivation of the classical 
RD equation. 

This high thermionic electron emission regime is ex¬ 
plained by the contribution of high energy electron 
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groups in the energy spectrum of the metal electrons. 
Therefore, Eq. m extends the predictions of the clas¬ 
sical Richardson-Dushman equation to higher metal tem¬ 
perature and permits to make us of k(T) to approximate 
the stationary energy distribution function of electrons 
in hot metals. 

Finally, regarding the applicability of these results, 
cathodes with higher thermionic electron emission cur¬ 
rents would led to improved characteristics in electron 
tube devices. Larger electron currents lead to greater 
amplification factors in traveling wave tubes or a better 


coupling with the microwave input signal in klystrons. 
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